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I. INTRODUCTION 

Quench a system from the homogeneous phase into a 
broken-symmetry one (think, e.g., to a ferromagnet or to 
a gas abruptly cooled below their critical temperature). 
The two phases have to separate and the process can be 
described via a local field u (local magnetization in ferro- 
magnets, density in liquid-vapor systems, concentration 
in alloys, etc.) on the physical space C K d . 

Very well known models pQ for the description of the 
field evolution are the Allen-Cahn and the Cahn-Hilliard 
equations. With suitable boundary conditions the former 
is an appropriate equation when the order parameter, i.e., 
the integral over fl of the field u, is not conserved, while 
the latter applies in the conserved case. On mathematical 
grounds, these equations can be thought as the gradient 
equation du/dt = —SF/Su, in a suitable Hilbert space 
[2^0] , for the Landau energy functional 



F[u] 



k 2 \\Vu\\ 2 



V(u) 



dx 



(1) 



with k > and the potential energy V a double well posi- 
tive regular function with degenerate zero value absolute 
minima in a, b, called phases of the system. More pre- 
cisely, if no constraint to the order parameter is imposed, 
it is possible to compute the gradient of the Landau func- 
tional in the Hilbert space L 2 (Vl) to get the Allen-Cahn 
equation [H equations (8) and (12)]. When the integral of 
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the field u is assumed to be constant throughout the evo- 
lution, computing the gradient in the L 2 (f2) results into a 
non-local evolution equation, while by using the Hilbert 
space ff _1 (f2) the Cahn-Hilliard equation is found [H 
equation (91) and the comment below equation (90)]. 

In this paper we discuss some properties of the station- 
ary profile connecting the two phases a and b, namely, the 
profile which is reached by the system at infinite time. 
More precisely, we explore the properties of the station- 
ary profile connecting the two phases when the double 
well potential energy is asymmetric. The equations for 
the interface profile are the Euler-Lagrange equations for 
the stationary points of the energy functional . These 
kind of problems are often referred to in the literature as 
gradient or diffuse interface problems [5]. 

We shall consider a one-field u system described by a 
potential energy V as above. More precisely we shall as- 
sume, throughout the paper, that the following condition 
is satisfied. 

Condition 1 The function V : K — » K is a positive 
C 2 (R) function with two single isolated local minima a 
and b (assume a < b) and such that V(a) = V(b) = 0, 
V"(a) > 0, and V"{b) > 0. 

Under this assumption, we shall approach the phase 
interface problem in dimension one, on a finite interval 
with Dirichlct boundary conditions fixing the values of 
the two phases at the extreme points of the interval to 
a and b, respectively. We prove that, for small fc, the 
position of the interface depends exclusively on the sec- 
ond derivatives of the potential energy V at the minima 
a and b. In particular, the statement is valid for the sta- 
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tionary solutions of the Allen-Cahn equation at imposed 
zero chemical potential and for those of the Cahn-Hilliard 
equation at imposed zero Lagrange multiplier (for mass 
conservation) . 

Note that even if V is not symmetric, the interface falls 
in the middle point of the space interval provided the two 
second derivatives are mutually equal. 

The interest to study the interface problem in the case 
of an asymmetric double-well potential energy is due to 
the fact that they appear naturally in different physi- 
cal situations. For example, in [3J HHH] we introduced a 
model describing a pressure guided transition between a 
fluid-rich and a fluid-poor phase and based on an asym- 
metric double-well potential energy. In this theory two 
order parameters m and e are introduced, having respec- 
tively the physical meaning of fluid density and solid 
strain. In [3J [7] the interface separating the two coex- 
isting phase has been studied and it has been shown that 
its localization properties, due to the asymmetry of the 
potential energy, are not trivial. 

Apart from the application mentioned above, due to 
its general character, our result is of interest for any sit- 
uation in which diffuse interfaces are relevant. In partic- 
ular when approximate computations or numerical sim- 
ulations are performed in the regime of small interface 
energy cost, our rigorous prediction could result valu- 
able in testing the soundness of the results. We remark 
that in this kind of problems numerical computations be- 
come particularly difficult when the interface energy cost 
is small (say k ~ I CP 2 with our parametrization) . 

The paper is organized as follows. In Section [IT] we 
state our general results. In Section III we discuss in 



detail their application to the porous media segregation 
problem. In Section |TV| we shall verify our result in some 
cases physically relevant by performing numerical compu- 
tations. In Section [VJ finally, we prove the results stated 
m Section HB 



II. RESULTS 

We approach the problem in dimension one. Recalling 
k is a positive constant, the Euler-Lagrange equations for 
the energy functional ([I]) reads as the Dirichlet boundary 
value problem 

k 2 u xx = V'(u) with u(0) = a and u{£) = b (2) 

for the field u(x), x g [0, £], for some I > 0. 




FIG. 1. Schematic description of the results in Theorems [T] 
and [2] The thick lines are the profiles at two finite values of 
k; the steepest one corresponds to the smallest value of k. 

From the physical point of view the problem ^ is that 
of finding the stationary profile connecting the two phases 
a and b on the finite space interval [0,£]. We remark 
that the above problem ^ is found when looking for 
the stationary solutions of the Allen-Cahn equation at 
imposed zero chemical potential or of the Cahn-Hilliard 
equation at imposed zero Lagrange multiplier (for mass 
conservation) with proper boundary conditions. 

Before stating our results, we note that by exploiting 
the one-dimensionality of the model a phase space anal- 
ysis [9] proves that the problem ^ has a unique solution 
implicitly given by the integral 

h ' ds (3) 



la ^2[E k + V{s)} 

where for any k > we have defined implicitly the energy 
level E k by the equation 



kds 



= £ 



(4) 



la yj2[E k + V( S ) 

namely, the integral in ([3]) with u = b and x — £. 

We note that the solution of ^ is a kink connecting on 
the interval [0, l\ the phase a to the phase b. At small k 
the energy cost associated with the gradient of u, see 
is small, so that the interface width is of order k and it 
is localized somewhere in the interval [0,£]. This remark 
is made rigorous in the next classical theorem, where we 
state that in the limit k — > the interface tends to a 
discontinuous kink connecting the two phases (see, also, 
figure [l]). 

Theorem 1 Assume Condition^is satisfied. Then, for 
any u\,U2 € (a, b) such that u\ < 112 

k ds 



lim 



k^oJ Ui ^2[E k + V(s) 



= 



(5) 
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FIG. 2. Schematic description of the heuristic interpretation of equation The thick line is the profile at finite k, whereas 
the thin discontinuous line is its k — > limit. 



The meaning of the theorem above is illustrated in 
figure [l] For any choice of a < u\ < U2 < b, by choosing 
k small enough the distance between the two points where 
the kink attains the values u\ and ui respectively can be 
made smaller than any fixed positive number. 

The theorem above suggests that for any k > 0, it 
is meaningful to define the interface position x\ as that 
(unique) point where the profile attains the value (a + 
b)/2, that is to say 



k di 



y/2[E k + V(s)} 



(6) 



Since for k — > the profile tends to a step function, dif- 
ferent definitions, equivalent in this limit, are possible 
for the interface position; e.g., the point where the pro- 
file equals the value of field corresponding to the local 
maximum of the potential energy. 

The following theorem gives the limiting behavior of 
the interface position for k —> 0. 

Theorem 2 Assume Condition^is satisfied. Then 

iim[vV>) 4 - (£ - 4)] = (7) 

fe— >o 

In words, the theorem states that in the limit k — > 
the distance of the interface from the boundary points 
is proportional to the inverse of the square root of the 
second derivative of the potential energy evaluated at 
the corresponding minimum of V. Note that from (JTJ) 
it follows immediately that the interface position x k for 
k — >• tends to 



(8) 



The limiting interface position Xg depends only on the 
second derivatives of the potential energy evaluated at 



the minima. In particular if V"(a) — V"(b) then x = 
1/2, that is to say, the interface is located at the middle 
point of the interval [0,£]. 

The proof (see Section [v| of Q is based on a direct 
evaluation of the integral Since is small for k — > 
the integrand diverges in a. The integral can be thus 
estimated by expanding in Taylor formula to the second 
order the potential energy V. 

Now, we give two heuristic arguments supporting the 



statement of the theorem. In Section IV on the other 



hand, we shall verify the result with some numerical com- 
putations in two cases of interest. 

Assuming k small, the kink is close to a step function. 
By phase space techniques it follows E k ~ 0. Fix rj > 
small and rj < 8 <C b — a. By ^ the distance between 
the points where the profile equals a + 77 and a + 5 is 



a+6 



k ds 



a+Tj 



^2[E k + V(s) 



a+S 



k ds 



VV^iu - a) 



where we have neglected E k and, noted that the integral 
is extend to a very small interval close to a, we expanded 
V in second order Taylor formula. Since rj is small, we 
have that the distance Lq +s from of the point where the 
value a + S is attained by the profile is approximatively 
given by L a + 5 a {k/y/W{j£)) log(<V*?)- 

By performing a similar computation in the neighbor- 
hood of the boundary point I and with obvious notation 
we get L b ~ s w (k/y/V"(b))\og(5/r)). The last two for- 
mulas allow to conjecture the validity of ([7]), see also 
figure [2j 

The second heuristic argument is based on the inter- 
pretation of ([2]) as the equation describing the motion of 
a particle of mass k 2 under the potential energy — V. In 
this language u is the position of the particle and x is 
the time. In the remaining part of this section, thus, we 
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shall address to u(x) as to the position of the particle at 
time x and, for this reason, derivatives with respect to x 
will be denoted by dots. 

In this context, solving ^ means looking for the mo- 
tion started at a with positive initial velocity vo such 
that at time £ the position b is reached. The conserva- 
tion of the total mechanical energy k 2 ii 2 /2 — V(u) implies 
that the motion will reach the point b with velocity vq- 
The interface position x\. is then the time at which the 
particle reaches the position (a + b)/2. 

We can describe such a motion by linearizing the equa- 
tion of motion in a neighborhood of the two unstable 
equilibrium points a and b. We get 



k 2 u = V"(a)(u-a) and k 2 u = V"{b){u - b) 



(b, —vq), respectively, and get 



u a (x) = a H sinh ( ^ 



for the motion started at a and 

u b (x) = b sinh ( — — x] 

^JWW) V k J 

for the motion that started at b. 

Recalling that x\ is the time at which the particle is 
at (a + b)/2, by the two equation above we get 



= arcsinhf 



b - a y/V"{a) 
2 kvo 



and 



k ^c:lnh( b - a " /V7W) 
^/Wib) V 2 kv Q 



respectively. We then solve the two equations with the 
initial conditions (u(0),u(0)) = (a, Vq) and (u(0),u(0)) = 



Recalling that arcsinh(.T) = log(x + \/x 2 + l), with 
simple algebra we get 



log 



y/V^ib -a) + ^V"{a){b-a) 2 +^k 2 v 2 



y/VWjib ~a) + ^V"{b){b-a) 2 + Ak 2 v 2 



(9) 



which suggests the validity of ([8]) since the right hand side 
tends to zero for k — >• 0. Note that vo does not diverge 
for k 0. 



III. PHASE INTERFACE IN POROUS MEDIA 

We discuss, now, in detail the application of the above 
results to the study of transitions from a fluid-poor to- 
wards a fluid-rich phase in porous media under consolida- 
tion. The model we study here is supposed to explain the 
behavior of the system in consolidation regime, namely, 
when an external pressure is applied on the solid compo- 
nent. The idea is that of explaining the existence of two 
phases, differing in liquid content, due to the modification 
in the solid structure caused by the external pressure. For 
this reason, as in the original Biot theory, in this model 
no chemical potential contribution is considered. 

In [31 [7] the interface separating the two coexisting 
phases has been studied numerically and it has been 
shown that its localization properties, due to the asym- 
metry of the potential energy, are not trivial. In fig- 



ures [3j[6] profiles for the solid strain e and the fluid den- 
sity m fields are shown. It is immediate to note that 
those kinks tend, for the interface energy cost tending to 
zero, to a not symmetric sharp interface between the two 
phases. 

More precisely, in this model three interface cost pa- 
rameters, fcxj ^2 j an d ^3, are introduced (see equation 
(111 below). These constants weights (e') 2 , e'm! 



(m') , respectively. The numerical simulations suggest 
that when ^1,^2,^3 — > with mutual ratios kept con- 
stant, the interface tends to a definite position not de- 
pending on these ratios. The e and the m profiles are 
shown in the cases k± = fea = kj, in figures [3J and g] and 
k% = 2/c 2 = 3&3 in figures [5] and [6j respectively. 

In Section [TT] we have proven results concerning the 
localization properties of a kink connecting two coexist- 
ing phases in a quite general one-field one-dimensional 
setup. As a straightforward application, we shall get the 
position of the limiting interface for the poromechanics 
model described above with the second gradient coeffi- 
cients such that k\k^ — k 2 = 0; in this particular case, 
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FIG. 3. Kink profiles for the strain field e for the po- 
tential energy ( |ll[ )-( fT3| ) with fei = k, k-z — k, k-j, = k, 
k = 10~\ 10~ 2 , 10~ 3 , a = 100, a = 1/2, and 6 = 1. 



FIG. 5. Kink profiles for the strain field e for the potential 
energy ([TT]> — ^13]) with ki = k, k 2 = k/2, k 3 = fe/3, k = 
10 _1 , 10~ 2 , 10" 3 , a = 100, a = 1/2, and 6=1. 
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FIG. 4. Kink profiles for the fluid density field m for the same 
parameters as in figure [3] 



FIG. 6. Kink profiles for the fluid density field m for the same 
parameters as in figure [5] 



indeed, the two-field model will reduce to a one-field one. 

For the complete description of the poromechanics 
setup we refer the readers to O [TU]. We recall 
that two fields m, e : [0,£\ — > K are introduced having 
the physical meaning of deviation from a reference value 
of the fluid density and solid strain, respectively. The 
stationary state of the system is described by the Euler- 
Lagrange equations, with suitable boundary conditions, 
for the variational problem associated with the Landau 
energy functional 



L[m,e}= [ [K(m',e') +^(m,s)]dx 
Jo 



(10) 



where K and ^ are, respectively, called the second and 
first gradient potential energy. 

The poromechanics model we study in this paper [5] 
is defined by choosing the second and the first gradient 
potential energies as follows. We let 

K(m',e r ) = i[fci(e') 2 + 2fc 2 eW + fc 3 (m') 2 ] (11) 



where fci, > and ^eK are such that kik^ — k\ > 
and the prime denotes the space derivative. Moreover we 
set 

*(m,£,p) = ^m 2 (3m 2 - 8bem + 6b 2 e 2 ) +^ B (m,e,p) 

(12) 
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where 



1 



1 



$ B (m,£;p) =pe + -e + -a(m - be) 



(13) 



is the Biot potential energy density [TT] , a > is the ratio 
between the fluid and the solid rigidity, b > is a cou- 
pling between the fluid and the solid component, p > 
is the external pressure, and a > is a material pa- 
rameter responsible for the showing up of the additional 
equilibrium. We remark that the condition k\k^ — k\ > 
ensures that the second gradient part K of the overall po- 
tential energy density is convex. Under this assumption 



there exists a minimizcr for the action functional ( 10 ) on 
a bounded domain [12 . 

In [SI |S] we have studied the minima of the potential 



() 



energy (12 1 and shown that they can be interpreted as 



the homogeneous phases of the system. We have proven 
that there exists a pressure p c , called critical pressure, 
such that for any p € [0,p c ) there exists a single phase 
(m s (p), e s (p)), called the standard phase, which is very 
similar to the usual solution of the Biot model. For 
p > p c a second phase (mf(p), £f(p)), richer in fluid with 
respect to the standard phase and hence called fluid- 
rich phase, appears. It has also been shown that for any 
p > the point (m s (p), £ s (f>)) is a local minimum of the 
two variable potential energy ty(m,e) with p fixed, while 
(mt(p), Sf(p)) is a local minimum for p > p c and a saddle 
point for p = p c . 

Moreover, in [7J it has been proven that there exists a 
unique value p co of the pressure, called coexistence pres- 
sure, such that the potential energy of the two phases 
is equal. More precisely, it has been proven that the 
equation ^ (m s (p) , e s (p)) = f(mf(p), £f{p)) has the sin- 
gle solution p co . We remark that it is not possible to 
provide explicit formulas for the phase field values and 
for the coexistence pressure. Indeed, these quantities do 
depend on the physical parameters of the model a, a, and 
b. In [7l H] we have proved (uniformly in these parame- 
ters) the scenario depicted above and we have provided 



a strategy for the numerical determination of all the in- 
teresting quantities. 

The behavior of the system at the coexistence pressure 
is particularly relevant; from now on we shall always con- 
sider p = p co and, for this reason, we shall drop p from 
the notation. When the external pressure is equal to p co , 
none of the two above phases is favored and we ask if 
profiles connecting one phase to the other exist [3] . In [7] 
this problem has been addressed on the space set R, while 
some results on a finite interval have been reported in [3] . 
There we have remarked that the localization properties 
of the kink profile are not trivial, as it has been fully 
explained above. 

The profiles connecting the two phases on a finite in- 
terval [0, £] at the coexistence pressure are given by the 
solutions (if any) of the Dirichlet boundary problem for 
the Euler-Lagrange equations associated with the func- 



tional (10), namely 



a* a ok _ Q and a* 

de dx de' dm 



d dK 
dx dm 1 



= 



(14) 



(see, for instance, [7J equation (4)]), with boundary con- 
ditions m(0) = m s , e(0) = e s , m(£) — rrif, and e(£) = £{. 
By writing explicitly the Euler-Lagrange equations we 
finally get the Dirichlet problem 



kie" + k 2 m" — — (2/3)abm 3 + ab 2 m 2 e + p + e — ab(m — be) 
k2e" + k^m" = am 3 — 2abm 2 e + ab 2 me 2 + a(m — be) 
m(0) = m s , e(0) = e s , rn{£) = irif, e(£) = £f 



(15) 



which can be approached numerically with the finite dif- 
ference method powered with the Newton-Raphson al- 
gorithm. The solution we find, in the cases of interest, 
are those depicted in the figures |3}|6| 



As noted above the behavior at small second gradi- 
ent coefficient fc 1; fc 2j an d fc 3 , of the stationary pro- 
files depicted in figures [3] - [6] is very peculiar: when 
k\, &2, ^3 — >■ with mutual ratios kept constant, the in- 
terface becomes sharp and its location tends to a definite 
position not depending on these ratios. Although the 
general results in Section In] do not apply to this case, 



we note that the solutions of the problem ( 15 ) behave in 



tune with those of the general one-field problem pi). 

At the moment, as already mentioned above, we can 
explain this fact only in the particular case kxk^ — fc| = 0, 
namely, fci/fe = ^2/^3- In this case, called degenerate 



case in [7J, the problem of finding a solution of the sta- 
tionary problem can be reduced to a one-field problem. 
Indeed, in such a case one performs the rotation of the 
Cartesian reference system 

m + Ae , — \m + e 



and 



77 := 



(16) 



VTTX 2 ' y/l + x 2 

in the plane m-e, where A := kxjk% = k 2 /k^, and defines 

U(t,r,) = V(m(Z,r ] ),e((,r l )) (17) 

Then one shows that the two fields m(x) and e(x) are 



solutions of the two equations ( 15 1 if and only if the cor- 



responding fields £(x) and rj(x) satisfy 
fe 3 (l + A 2 )£"=- (i.;/) and 



dU (t \ A dU (C ^ 
(£,?/) and ^-(£,?7) 



(18) 



9£ drj 

The root locus of the constraint curve dU(£, r])/dr] = 
is made of a certain number of maximal components 
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such that each of them is the graph of a function £ £ 
K — > 7j(£) G K; for each of them the first between the 



two equations ( 18 ) becomes a one-field one-dimensional 
problem. 

Since the function U has been obtained by rotating 
the coordinate axes, then at the coexistence pressure it 
has the two absolute minimum points (£ s ,?7s) and (£f, r?f) 
corresponding, respectively, to the standard and to the 
fluid-rich phases. Since (m s ,e s ) and (mf, £f) satisfy the 
equations v E , m (m, e) = and ^/ E (m, e) = 0, we have that 
the two points (£ s , r] s ) and (£f, rjf) are solutions of the con- 
straint equation dU(£,rj)/dr) = and hence they belong 
to the constraint curve. 

In [7] we have seen that there exist values of the second 
gradient parameters k\, Jc2, and k% such that the two 
points above fall on the same maximal component of the 
constraint equation. Since, in this case, the function U 
has two isolated absolute minima which, by hypothesis, 
belong to the same maximal component of the constraint 
curve, we have that the function £/(£, ??(£)) °f £ has two 
absolute isolated minima in £ s and £f. 

In this particular case, finally, the theory developed in 
Section III] can be applied to the poromcchanics problem 



(15) and the results discussed in the introduction find a 



complete explanation. The explanation is not only qual- 
itative, but also quantitative, indeed, for the first gradi- 
ent parameters chosen as in the figures [3] and |4j namely, 
a = 100, a = 1/2, and b = 1, by computing the second 
derivatives of the one-field potential energy U (£, at 
the phases £ s and £f and by using Q we get 0.6164 for 
the position of the interface on the interval [0, 1]. For the 
second gradient parameter choice in figures [3] and [4] the 



degeneracy condition kik^ 



is satisfied. So that 



in this case the theory developed in Section [TT| applies 
and the match with the numerical result is striking. 



IV. TWO NUMERICAL EXAMPLES 

In this section we check numerically |7| for two very 
well known potentials. The first example, see (19), is 



an asymmetric potential with mutually different second 
derivatives at the minima; the numerical computation 
will confirm that, for k small, the interface is not located 
at the middle point of the space interval and that its posi- 
tion is given by (fsb - The second example, see (20), is, on 



the other hand, an asymmetric potential with mutually 
equal second derivatives at the minima; the numerical 
computation will confirm that, for k small, the interface 



is located at the middle point of the space interval. 

We shall solve the equation ^ with £ = 1 via the 
following procedure. Fix k and use Q to compute the 
energy level E^; then use ^ to compute the profile. 

We consider first the following asymmetric double-well 
potential [T31IT4] 



^i(u) = 



1 / 1 — exp[&i(u — 1)] 
1 - exp[6i(u - 1)] 
l-exp[-6 2 (u+l)] 



x 



1 - exp[-& 2 (wo + 1)] 



(19) 



where uj £ R, 61,62 > 0, and u £ (— 1,+1). The 
two minima are at ±1 and uq controls the height of 
the internal barrier. The kinks corresponding to k = 
1.0,0.5,0.25,0.1 (for u = 1, b x = 1, b 2 = 5, u = 0.5) 
are shown in figure [7J By applying ^ one gets Xq = 
[1 + e 2 + e 4 + e 6 + e 8 ]/[l + e 2 + e 4 + e 6 + 6e 8 ] » 0.187846, 
so that the numerical computation is in good agreement 
with the theoretical prediction (j8|). 

Another interesting example is the asymmetric poten- 
tial which is considered the benchmark potential [T5HT7] 
of the current literature on rocked ratchets, namely, 



V 2 {u) 



2tt 



(u-u ) 



47T 



(u-Uq) 



-R (20) 



with a £ R, a- 1 = (a/(8 7 r))(V3/2) 1 / 2 (3 + v^), u = 
(a/(27r))cos- 1 ((-l + V3)/2), and R = -2ir/a. The not 



essential additive constant R in ( 20 ) has been introduced 



in order to ensure that the potential is equal to zero at 
the minima na with n £ Z. The second derivative of 
V2 at the minima na are all equal since the function is 
periodic with period a. Some algebra yields F 2 "( M ) = 
16(— l + v / 3)7r 3 /a 3 , so that these derivatives are not zero. 
The kinks corresponding to k = 3.0, 1.5, 1.0, 0.5, 0.45 (for 
a = 1) are shown in figure [8j 

We stress that, in this case, although the potential en- 
ergy is asymmetric with respect to the phase exchange, 
the interface tends to be located in the middle point of 
the space interval, since the second derivatives of the po- 
tential energy at the minima are mutually equal. 

Note that our result applies even if deformations of the 



standard rocket ratchet potential ( 20 ) are considered (sec 
[18] ) . For instance the positive multiplicative factor 1/4 



in the second term in ( 20 ) can be chosen differently. Our 



result is still valid, provided no other minima is found 
between and a = 1, since the second derivatives of the 
potential at the minima are equal and do not vanish. 

It is notable to remark the following. In the above 
sections we have studied the localization properties of 
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FIG. 7. Graph of the kink for the potential (19 1 for cjo = 1, 
bi = 1, 62 = 5, and no = 0.5. The profiles correspond to 
k — 1.0,0.5,0.25,0.1 in increasing order of steepness. 



the interfaces connecting two phases under an asymmet- 
ric double-well potential energy with degenerate minima 
and we have stated (see Section [V] for the proof) that 
the position of the interface depends exclusively on the 
second derivatives of the potential energies at the two 
phases. 

Although we cannot say, at the moment, to which ex- 
tent this analogy between the two results can be pushed 
forward, we note that our result is similar in spirit to 
those related to the direction of motion of a dc soliton 
current in ac drive rocket ratchet. Which is, a priori, a 
completely different problem with respect to the phase 
interface problem ^ that we have studied. 

We refer in particular to [13j [14] . There the authors 



FIG. 8. Graph of the kink for the potential (20 1 for a = 1. The 



considered the potential ( 19 1 and proved the following: a 



string soliton profile is pushed by "symmetric" thermal 
fluctuations sidewise towards the boundary point where 
the boundary condition has been chosen equal to the 
phase corresponding to the narrowest valley (biggest sec- 
ond derivative) of the double-well potential. In particu- 
lar in |14j it is shown that this tendency to move towards 
the narrowest valley is controlled by the second deriva- 
tives of the potential energy evaluated at the phases. 

V. PROOF OF RESULTS 

The proof of the results in Section [IT] is based on 
the classical Weierstrass qualitative analysis of the one- 
dimensional mechanical problem equivalent to ^ [S]. 

Before starting the proof of the results, we note that E k 
is an increasing positive function of k and E k — > for k — > 
0. This is a quite intuitive result since taking the limit 
k — > is the same as considering a zero mass limit in the 



profiles correspond to k = 3.0, 1.5, 1.0,0.5,0.45 in increasing 
order of steepness. The two profiles corresponding to k = 0.5 
and k = 0.45 are almost coincident in the picture. 

equivalent mechanical system. More precisely, remark 
that 



a y/2[E + V(s)] 



is an decreasing function of E, defined for E > 0, and 
such that its limit for E — > tends to infinity. Since by 



ds 



a y/2[E~k + V(s)} k 

it follows that for k — > the energy level Ek tends to 0. 

The Theorem [T] is a classical result in interface theory. 
For the sake of completeness we report the proof in our 
case. 

Proof of Theorem^ Since V is a strictly positive func- 
tion in [m, U2], by the dominated convergence theorem it 
follows immediately that 

<" U2 ds r U2 ds 



lim 



k^oJ Ui ^2[E k + 7(8)] J U1 y/2V{s) 
yielding |5]). 



< 00 



□ 



The proof of Theorem [2] based on the heuristic argu- 
ments discussed in Section [TTJ uses some ideas developed 
in [ini Appendix A]. 

Proof of Theorem [^[ Recall the definition (Jf)]) of x\ . 
By simple algebra we have that 



where 



S a ., 



ds 



(21) 



y/2[E k + V"(a)(8-a) 2 /2] 
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and 



ary condition b has been fixed. Noted that 



R a ,k — 



ds 



^2[E k + V{sj] 



By direct integration we get 



S a ,k — 



arcsinh 



'V'(a) b- a 



2E k 



(22) 



The idea behind the proof is the following: once S a , k 
has been computed explicitly, we are able to manage the 
divergence of the integral with E k for k — > 0. On the 
other hand we shall simply estimate R a ,k and prove that 
it is finite for any k > 0. 

We can perform the analogous computation in the 
neighborhood of the boundary point I where the bound- 



we get 



where 



and 



k ds 



±> y/2[E k + V(s)} 



Sb,k = 



£- x 1 k = kS b .k + kR b , k 

b ds 

±i y/2[E k + V"(b){b- s) 2 /2] 

1 



arcsinh 



! V"(b) b-a 
2E k 



R 



b.k 



ds 



Sb,k 



(23) 



(24) 



ts y/2[E k + V(s)} 
Recalling that arcsinh(x) = log(a; + v / x 2 ~+l) , by equa- 



tions (21)— ( 24) we get 



JV^)x k -JW(Vj{t-x k ) = k 



log 



y/V*(a)(b-a) + ^V"{a){b-ay + %E k ATTTmp 

/T77777T/, n , H r „nsn + V V W R a,k ~ V V ( b ) R b,k 



y/V*(b)(b - a) + y/V"(b)(b - a) 2 + 8E k 



(25) 



which is the analogous of the equation (|9| that we ob- 
tained in the heuristic discussion in Section HTl 



By using the explicit expression of S a ,k and Sb, k we 
have been able to cancel the divergence in E k in the first 
part of the right-hand side of (25). Recall equation Q; 



since E k — > for k — > 0, the theorem will follow from 



( 25 ) once we shall have proven that 



sup|i? a .fc| < oo and sup|i?b,fc| < co (26) 

fc>0 fc>0 

We prove the first of the two bounds above; the second 
one can be proven similarly. With some algebra we get 



R a ,, 



y/2E k + V"(a)(s - a) 2 - y/2[E k - V(s)} ^ 
^/2E k + V"{a){s - a)^2[E k - V(s)} 

[-2V(s) + V"{a)(s-a) 2 }ds 

y/2[E k + V(s)][2E k + V"(a)(s - a) 2 ] + 2[E k + V(s)\ ^2E k + V"(a){s - af 



where in the last step we have multiplied numerator and Since, as we remarked at the beginning of this section, 
denominator times the sum of the two square root terms. E k > for any k > and E k — > for k — > 0, we have 

that 
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SUp \R a ,k\ < 
k>0 



\-2V(s) + V"{a){s-a) 2 \ds 
2V(s)^V"(a)(s -a) + y/2V(s)V"(a)(s - a) 2 



Moreover the integrand is continuous in (a, (a + 6)/2] and 
its limit for s -> a is \V"'{a)\/[6V" (a)y/V"{a)]; therefore, 
we have that the integral on the right-hand side of the 
inequality above is finite. We thus get the first of the two 



bounds (26) 



As noted below (25) this is sufficient to complete the 



proof of the theorem. 



□ 
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